3.3.51 \(\int \frac {1}{\sqrt {a+\frac {b}{x}} (c+\frac {d}{x})^3} \, dx\) [251]
Optimal. Leaf size=250 \[ \frac {d (2 b c-3 a d) \sqrt {a+\frac {b}{x}}}{2 a c^2 (b c-a d) \left (c+\frac {d}{x}\right )^2}+\frac {d (b c-4 a d) (4 b c-3 a d) \sqrt {a+\frac {b}{x}}}{4 a c^3 (b c-a d)^2 \left (c+\frac {d}{x}\right )}+\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}-\frac {d^{3/2} \left (35 b^2 c^2-56 a b c d+24 a^2 d^2\right ) \tan ^{-1}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b c-a d}}\right )}{4 c^4 (b c-a d)^{5/2}}-\frac {(b c+6 a d) \tanh ^{-1}\left (\frac {\sqrt {a+\frac {b}{x}}}{\sqrt {a}}\right )}{a^{3/2} c^4} \]
[Out]
-1/4*d^(3/2)*(24*a^2*d^2-56*a*b*c*d+35*b^2*c^2)*arctan(d^(1/2)*(a+b/x)^(1/2)/(-a*d+b*c)^(1/2))/c^4/(-a*d+b*c)^
(5/2)-(6*a*d+b*c)*arctanh((a+b/x)^(1/2)/a^(1/2))/a^(3/2)/c^4+1/2*d*(-3*a*d+2*b*c)*(a+b/x)^(1/2)/a/c^2/(-a*d+b*
c)/(c+d/x)^2+1/4*d*(-4*a*d+b*c)*(-3*a*d+4*b*c)*(a+b/x)^(1/2)/a/c^3/(-a*d+b*c)^2/(c+d/x)+x*(a+b/x)^(1/2)/a/c/(c
+d/x)^2
________________________________________________________________________________________
Rubi [A]
time = 0.25, antiderivative size = 250, normalized size of antiderivative = 1.00, number of steps
used = 9, number of rules used = 7, integrand size = 21, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.333, Rules used = {382, 105, 156,
162, 65, 214, 211} \begin {gather*} -\frac {(6 a d+b c) \tanh ^{-1}\left (\frac {\sqrt {a+\frac {b}{x}}}{\sqrt {a}}\right )}{a^{3/2} c^4}-\frac {d^{3/2} \left (24 a^2 d^2-56 a b c d+35 b^2 c^2\right ) \text {ArcTan}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b c-a d}}\right )}{4 c^4 (b c-a d)^{5/2}}+\frac {d \sqrt {a+\frac {b}{x}} (b c-4 a d) (4 b c-3 a d)}{4 a c^3 \left (c+\frac {d}{x}\right ) (b c-a d)^2}+\frac {d \sqrt {a+\frac {b}{x}} (2 b c-3 a d)}{2 a c^2 \left (c+\frac {d}{x}\right )^2 (b c-a d)}+\frac {x \sqrt {a+\frac {b}{x}}}{a c \left (c+\frac {d}{x}\right )^2} \end {gather*}
Antiderivative was successfully verified.
[In]
Int[1/(Sqrt[a + b/x]*(c + d/x)^3),x]
[Out]
(d*(2*b*c - 3*a*d)*Sqrt[a + b/x])/(2*a*c^2*(b*c - a*d)*(c + d/x)^2) + (d*(b*c - 4*a*d)*(4*b*c - 3*a*d)*Sqrt[a
+ b/x])/(4*a*c^3*(b*c - a*d)^2*(c + d/x)) + (Sqrt[a + b/x]*x)/(a*c*(c + d/x)^2) - (d^(3/2)*(35*b^2*c^2 - 56*a*
b*c*d + 24*a^2*d^2)*ArcTan[(Sqrt[d]*Sqrt[a + b/x])/Sqrt[b*c - a*d]])/(4*c^4*(b*c - a*d)^(5/2)) - ((b*c + 6*a*d
)*ArcTanh[Sqrt[a + b/x]/Sqrt[a]])/(a^(3/2)*c^4)
Rule 65
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_), x_Symbol] :> With[{p = Denominator[m]}, Dist[p/b, Sub
st[Int[x^(p*(m + 1) - 1)*(c - a*(d/b) + d*(x^p/b))^n, x], x, (a + b*x)^(1/p)], x]] /; FreeQ[{a, b, c, d}, x] &
& NeQ[b*c - a*d, 0] && LtQ[-1, m, 0] && LeQ[-1, n, 0] && LeQ[Denominator[n], Denominator[m]] && IntLinearQ[a,
b, c, d, m, n, x]
Rule 105
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_.), x_Symbol] :> Simp[b*(a +
b*x)^(m + 1)*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/((m + 1)*(b*c - a*d)*(b*e - a*f))), x] + Dist[1/((m + 1)*(b*
c - a*d)*(b*e - a*f)), Int[(a + b*x)^(m + 1)*(c + d*x)^n*(e + f*x)^p*Simp[a*d*f*(m + 1) - b*(d*e*(m + n + 2) +
c*f*(m + p + 2)) - b*d*f*(m + n + p + 3)*x, x], x], x] /; FreeQ[{a, b, c, d, e, f, n, p}, x] && ILtQ[m, -1] &
& (IntegerQ[n] || IntegersQ[2*n, 2*p] || ILtQ[m + n + p + 3, 0])
Rule 156
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_))^(p_)*((g_.) + (h_.)*(x_)), x_Symb
ol] :> Simp[(b*g - a*h)*(a + b*x)^(m + 1)*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/((m + 1)*(b*c - a*d)*(b*e - a*f
))), x] + Dist[1/((m + 1)*(b*c - a*d)*(b*e - a*f)), Int[(a + b*x)^(m + 1)*(c + d*x)^n*(e + f*x)^p*Simp[(a*d*f*
g - b*(d*e + c*f)*g + b*c*e*h)*(m + 1) - (b*g - a*h)*(d*e*(n + 1) + c*f*(p + 1)) - d*f*(b*g - a*h)*(m + n + p
+ 3)*x, x], x], x] /; FreeQ[{a, b, c, d, e, f, g, h, n, p}, x] && ILtQ[m, -1]
Rule 162
Int[(((e_.) + (f_.)*(x_))^(p_)*((g_.) + (h_.)*(x_)))/(((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))), x_Symbol] :>
Dist[(b*g - a*h)/(b*c - a*d), Int[(e + f*x)^p/(a + b*x), x], x] - Dist[(d*g - c*h)/(b*c - a*d), Int[(e + f*x)
^p/(c + d*x), x], x] /; FreeQ[{a, b, c, d, e, f, g, h}, x]
Rule 211
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]
Rule 214
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x/Rt[-a/b, 2]], x] /; FreeQ[{a, b},
x] && NegQ[a/b]
Rule 382
Int[((a_) + (b_.)*(x_)^(n_))^(p_.)*((c_) + (d_.)*(x_)^(n_))^(q_.), x_Symbol] :> -Subst[Int[(a + b/x^n)^p*((c +
d/x^n)^q/x^2), x], x, 1/x] /; FreeQ[{a, b, c, d, p, q}, x] && NeQ[b*c - a*d, 0] && ILtQ[n, 0]
Rubi steps
\begin {align*} \int \frac {1}{\sqrt {a+\frac {b}{x}} \left (c+\frac {d}{x}\right )^3} \, dx &=-\text {Subst}\left (\int \frac {1}{x^2 \sqrt {a+b x} (c+d x)^3} \, dx,x,\frac {1}{x}\right )\\ &=\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}+\frac {\text {Subst}\left (\int \frac {\frac {1}{2} (b c+6 a d)+\frac {5 b d x}{2}}{x \sqrt {a+b x} (c+d x)^3} \, dx,x,\frac {1}{x}\right )}{a c}\\ &=\frac {d (2 b c-3 a d) \sqrt {a+\frac {b}{x}}}{2 a c^2 (b c-a d) \left (c+\frac {d}{x}\right )^2}+\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}-\frac {\text {Subst}\left (\int \frac {-(b c-a d) (b c+6 a d)-\frac {3}{2} b d (2 b c-3 a d) x}{x \sqrt {a+b x} (c+d x)^2} \, dx,x,\frac {1}{x}\right )}{2 a c^2 (b c-a d)}\\ &=\frac {d (2 b c-3 a d) \sqrt {a+\frac {b}{x}}}{2 a c^2 (b c-a d) \left (c+\frac {d}{x}\right )^2}+\frac {d (b c-4 a d) (4 b c-3 a d) \sqrt {a+\frac {b}{x}}}{4 a c^3 (b c-a d)^2 \left (c+\frac {d}{x}\right )}+\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}+\frac {\text {Subst}\left (\int \frac {(b c-a d)^2 (b c+6 a d)+\frac {1}{4} b d (b c-4 a d) (4 b c-3 a d) x}{x \sqrt {a+b x} (c+d x)} \, dx,x,\frac {1}{x}\right )}{2 a c^3 (b c-a d)^2}\\ &=\frac {d (2 b c-3 a d) \sqrt {a+\frac {b}{x}}}{2 a c^2 (b c-a d) \left (c+\frac {d}{x}\right )^2}+\frac {d (b c-4 a d) (4 b c-3 a d) \sqrt {a+\frac {b}{x}}}{4 a c^3 (b c-a d)^2 \left (c+\frac {d}{x}\right )}+\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}+\frac {(b c+6 a d) \text {Subst}\left (\int \frac {1}{x \sqrt {a+b x}} \, dx,x,\frac {1}{x}\right )}{2 a c^4}-\frac {\left (d^2 \left (35 b^2 c^2-56 a b c d+24 a^2 d^2\right )\right ) \text {Subst}\left (\int \frac {1}{\sqrt {a+b x} (c+d x)} \, dx,x,\frac {1}{x}\right )}{8 c^4 (b c-a d)^2}\\ &=\frac {d (2 b c-3 a d) \sqrt {a+\frac {b}{x}}}{2 a c^2 (b c-a d) \left (c+\frac {d}{x}\right )^2}+\frac {d (b c-4 a d) (4 b c-3 a d) \sqrt {a+\frac {b}{x}}}{4 a c^3 (b c-a d)^2 \left (c+\frac {d}{x}\right )}+\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}+\frac {(b c+6 a d) \text {Subst}\left (\int \frac {1}{-\frac {a}{b}+\frac {x^2}{b}} \, dx,x,\sqrt {a+\frac {b}{x}}\right )}{a b c^4}-\frac {\left (d^2 \left (35 b^2 c^2-56 a b c d+24 a^2 d^2\right )\right ) \text {Subst}\left (\int \frac {1}{c-\frac {a d}{b}+\frac {d x^2}{b}} \, dx,x,\sqrt {a+\frac {b}{x}}\right )}{4 b c^4 (b c-a d)^2}\\ &=\frac {d (2 b c-3 a d) \sqrt {a+\frac {b}{x}}}{2 a c^2 (b c-a d) \left (c+\frac {d}{x}\right )^2}+\frac {d (b c-4 a d) (4 b c-3 a d) \sqrt {a+\frac {b}{x}}}{4 a c^3 (b c-a d)^2 \left (c+\frac {d}{x}\right )}+\frac {\sqrt {a+\frac {b}{x}} x}{a c \left (c+\frac {d}{x}\right )^2}-\frac {d^{3/2} \left (35 b^2 c^2-56 a b c d+24 a^2 d^2\right ) \tan ^{-1}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b c-a d}}\right )}{4 c^4 (b c-a d)^{5/2}}-\frac {(b c+6 a d) \tanh ^{-1}\left (\frac {\sqrt {a+\frac {b}{x}}}{\sqrt {a}}\right )}{a^{3/2} c^4}\\ \end {align*}
________________________________________________________________________________________
Mathematica [A]
time = 1.25, size = 215, normalized size = 0.86 \begin {gather*} \frac {\frac {c \sqrt {a+\frac {b}{x}} x \left (4 b^2 c^2 (d+c x)^2+2 a^2 d^2 \left (6 d^2+9 c d x+2 c^2 x^2\right )-a b c d \left (19 d^2+29 c d x+8 c^2 x^2\right )\right )}{a (b c-a d)^2 (d+c x)^2}-\frac {d^{3/2} \left (35 b^2 c^2-56 a b c d+24 a^2 d^2\right ) \tan ^{-1}\left (\frac {\sqrt {d} \sqrt {a+\frac {b}{x}}}{\sqrt {b c-a d}}\right )}{(b c-a d)^{5/2}}-\frac {4 (b c+6 a d) \tanh ^{-1}\left (\frac {\sqrt {a+\frac {b}{x}}}{\sqrt {a}}\right )}{a^{3/2}}}{4 c^4} \end {gather*}
Antiderivative was successfully verified.
[In]
Integrate[1/(Sqrt[a + b/x]*(c + d/x)^3),x]
[Out]
((c*Sqrt[a + b/x]*x*(4*b^2*c^2*(d + c*x)^2 + 2*a^2*d^2*(6*d^2 + 9*c*d*x + 2*c^2*x^2) - a*b*c*d*(19*d^2 + 29*c*
d*x + 8*c^2*x^2)))/(a*(b*c - a*d)^2*(d + c*x)^2) - (d^(3/2)*(35*b^2*c^2 - 56*a*b*c*d + 24*a^2*d^2)*ArcTan[(Sqr
t[d]*Sqrt[a + b/x])/Sqrt[b*c - a*d]])/(b*c - a*d)^(5/2) - (4*(b*c + 6*a*d)*ArcTanh[Sqrt[a + b/x]/Sqrt[a]])/a^(
3/2))/(4*c^4)
________________________________________________________________________________________
Maple [B] Leaf count of result is larger than twice the leaf count of optimal. \(2268\) vs.
\(2(222)=444\).
time = 0.10, size = 2269, normalized size = 9.08
| | |
method |
result |
size |
| | |
risch |
\(\text {Expression too large to display}\) |
\(1288\) |
default |
\(\text {Expression too large to display}\) |
\(2269\) |
| | |
|
|
|
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(1/(c+d/x)^3/(a+1/x*b)^(1/2),x,method=_RETURNVERBOSE)
[Out]
1/8*((a*x+b)/x)^(1/2)*x*(-91*a^(7/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x
+d))*b^2*c^2*d^5+35*a^(5/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b^3*
c^3*d^4+24*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^2*(d*(a*d-b*c)/c^2)^(1/2)*b^3*c^5*d^2*x+46*
(x*(a*x+b))^(1/2)*a^(5/2)*(d*(a*d-b*c)/c^2)^(1/2)*b^2*c^6*d*x^2+68*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b
)/a^(1/2))*a^4*(d*(a*d-b*c)/c^2)^(1/2)*b*c^4*d^3*x^2-60*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*
a^3*(d*(a*d-b*c)/c^2)^(1/2)*b^2*c^5*d^2*x^2+12*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^2*(d*(a
*d-b*c)/c^2)^(1/2)*b^3*c^6*d*x^2-102*(x*(a*x+b))^(1/2)*a^(7/2)*(d*(a*d-b*c)/c^2)^(1/2)*b*c^4*d^3*x+92*(x*(a*x+
b))^(1/2)*a^(5/2)*(d*(a*d-b*c)/c^2)^(1/2)*b^2*c^5*d^2*x+136*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/
2))*a^4*(d*(a*d-b*c)/c^2)^(1/2)*b*c^3*d^4*x-120*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^3*(d*(
a*d-b*c)/c^2)^(1/2)*b^2*c^4*d^3*x+22*(x*(a*x+b))^(1/2)*a^(7/2)*(d*(a*d-b*c)/c^2)^(1/2)*b*c^6*d*x^3-22*(x*(a*x+
b))^(3/2)*a^(5/2)*(d*(a*d-b*c)/c^2)^(1/2)*b*c^6*d*x-18*(x*(a*x+b))^(1/2)*a^(7/2)*(d*(a*d-b*c)/c^2)^(1/2)*b*c^5
*d^2*x^2-12*(x*(a*x+b))^(1/2)*a^(9/2)*(d*(a*d-b*c)/c^2)^(1/2)*c^5*d^2*x^3+12*(x*(a*x+b))^(3/2)*a^(7/2)*(d*(a*d
-b*c)/c^2)^(1/2)*c^5*d^2*x-24*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^5*(d*(a*d-b*c)/c^2)^(1/2
)*c^3*d^4*x^2+80*a^(9/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b*c^3*d
^4*x^2-91*a^(7/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b^2*c^4*d^3*x^
2+35*a^(5/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b^3*c^5*d^2*x^2-18*
(x*(a*x+b))^(3/2)*a^(5/2)*(d*(a*d-b*c)/c^2)^(1/2)*b*c^5*d^2+36*(x*(a*x+b))^(1/2)*a^(9/2)*(d*(a*d-b*c)/c^2)^(1/
2)*c^3*d^4*x-16*(x*(a*x+b))^(1/2)*a^(3/2)*(d*(a*d-b*c)/c^2)^(1/2)*b^3*c^6*d*x-48*ln(1/2*(2*(x*(a*x+b))^(1/2)*a
^(1/2)+2*a*x+b)/a^(1/2))*a^5*(d*(a*d-b*c)/c^2)^(1/2)*c^2*d^5*x+8*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/
a^(1/2))*a*(d*(a*d-b*c)/c^2)^(1/2)*b^4*c^6*d*x+160*a^(9/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2
*a*d*x+b*c*x-b*d)/(c*x+d))*b*c^2*d^5*x-8*(x*(a*x+b))^(1/2)*a^(3/2)*(d*(a*d-b*c)/c^2)^(1/2)*b^3*c^7*x^2-182*a^(
7/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b^2*c^3*d^4*x+70*a^(5/2)*ln
((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b^3*c^4*d^3*x-62*(x*(a*x+b))^(1/2)
*a^(7/2)*(d*(a*d-b*c)/c^2)^(1/2)*b*c^3*d^4+46*(x*(a*x+b))^(1/2)*a^(5/2)*(d*(a*d-b*c)/c^2)^(1/2)*b^2*c^4*d^3+68
*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^4*(d*(a*d-b*c)/c^2)^(1/2)*b*c^2*d^5-60*ln(1/2*(2*(x*(
a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^3*(d*(a*d-b*c)/c^2)^(1/2)*b^2*c^3*d^4+12*ln(1/2*(2*(x*(a*x+b))^(1/2)
*a^(1/2)+2*a*x+b)/a^(1/2))*a^2*(d*(a*d-b*c)/c^2)^(1/2)*b^3*c^4*d^3+4*ln(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x
+b)/a^(1/2))*a*(d*(a*d-b*c)/c^2)^(1/2)*b^4*c^7*x^2-24*a^(11/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)
*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*c^2*d^5*x^2+8*(x*(a*x+b))^(3/2)*a^(7/2)*(d*(a*d-b*c)/c^2)^(1/2)*c^4*d^3-48*a^(1
1/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*c*d^6*x+24*(x*(a*x+b))^(1/2
)*a^(9/2)*(d*(a*d-b*c)/c^2)^(1/2)*c^2*d^5-8*(x*(a*x+b))^(1/2)*a^(3/2)*(d*(a*d-b*c)/c^2)^(1/2)*b^3*c^5*d^2-24*l
n(1/2*(2*(x*(a*x+b))^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a^5*(d*(a*d-b*c)/c^2)^(1/2)*c*d^6+4*ln(1/2*(2*(x*(a*x+b))
^(1/2)*a^(1/2)+2*a*x+b)/a^(1/2))*a*(d*(a*d-b*c)/c^2)^(1/2)*b^4*c^5*d^2+80*a^(9/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(
a*d-b*c)/c^2)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*b*c*d^6-24*a^(11/2)*ln((2*(x*(a*x+b))^(1/2)*(d*(a*d-b*c)/c^2
)^(1/2)*c-2*a*d*x+b*c*x-b*d)/(c*x+d))*d^7)/c^5/(x*(a*x+b))^(1/2)/(a*d-b*c)^3/(c*x+d)^2/a^(5/2)/(d*(a*d-b*c)/c^
2)^(1/2)
________________________________________________________________________________________
Maxima [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \text {Failed to integrate} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c+d/x)^3/(a+b/x)^(1/2),x, algorithm="maxima")
[Out]
integrate(1/(sqrt(a + b/x)*(c + d/x)^3), x)
________________________________________________________________________________________
Fricas [B] Leaf count of result is larger than twice the leaf count of optimal. 568 vs.
\(2 (222) = 444\).
time = 4.36, size = 2307, normalized size = 9.23 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c+d/x)^3/(a+b/x)^(1/2),x, algorithm="fricas")
[Out]
[1/8*(4*(b^3*c^3*d^2 + 4*a*b^2*c^2*d^3 - 11*a^2*b*c*d^4 + 6*a^3*d^5 + (b^3*c^5 + 4*a*b^2*c^4*d - 11*a^2*b*c^3*
d^2 + 6*a^3*c^2*d^3)*x^2 + 2*(b^3*c^4*d + 4*a*b^2*c^3*d^2 - 11*a^2*b*c^2*d^3 + 6*a^3*c*d^4)*x)*sqrt(a)*log(2*a
*x - 2*sqrt(a)*x*sqrt((a*x + b)/x) + b) + (35*a^2*b^2*c^2*d^3 - 56*a^3*b*c*d^4 + 24*a^4*d^5 + (35*a^2*b^2*c^4*
d - 56*a^3*b*c^3*d^2 + 24*a^4*c^2*d^3)*x^2 + 2*(35*a^2*b^2*c^3*d^2 - 56*a^3*b*c^2*d^3 + 24*a^4*c*d^4)*x)*sqrt(
-d/(b*c - a*d))*log(-(2*(b*c - a*d)*x*sqrt(-d/(b*c - a*d))*sqrt((a*x + b)/x) - b*d + (b*c - 2*a*d)*x)/(c*x + d
)) + 2*(4*(a*b^2*c^5 - 2*a^2*b*c^4*d + a^3*c^3*d^2)*x^3 + (8*a*b^2*c^4*d - 29*a^2*b*c^3*d^2 + 18*a^3*c^2*d^3)*
x^2 + (4*a*b^2*c^3*d^2 - 19*a^2*b*c^2*d^3 + 12*a^3*c*d^4)*x)*sqrt((a*x + b)/x))/(a^2*b^2*c^6*d^2 - 2*a^3*b*c^5
*d^3 + a^4*c^4*d^4 + (a^2*b^2*c^8 - 2*a^3*b*c^7*d + a^4*c^6*d^2)*x^2 + 2*(a^2*b^2*c^7*d - 2*a^3*b*c^6*d^2 + a^
4*c^5*d^3)*x), 1/8*(8*(b^3*c^3*d^2 + 4*a*b^2*c^2*d^3 - 11*a^2*b*c*d^4 + 6*a^3*d^5 + (b^3*c^5 + 4*a*b^2*c^4*d -
11*a^2*b*c^3*d^2 + 6*a^3*c^2*d^3)*x^2 + 2*(b^3*c^4*d + 4*a*b^2*c^3*d^2 - 11*a^2*b*c^2*d^3 + 6*a^3*c*d^4)*x)*s
qrt(-a)*arctan(sqrt(-a)*sqrt((a*x + b)/x)/a) + (35*a^2*b^2*c^2*d^3 - 56*a^3*b*c*d^4 + 24*a^4*d^5 + (35*a^2*b^2
*c^4*d - 56*a^3*b*c^3*d^2 + 24*a^4*c^2*d^3)*x^2 + 2*(35*a^2*b^2*c^3*d^2 - 56*a^3*b*c^2*d^3 + 24*a^4*c*d^4)*x)*
sqrt(-d/(b*c - a*d))*log(-(2*(b*c - a*d)*x*sqrt(-d/(b*c - a*d))*sqrt((a*x + b)/x) - b*d + (b*c - 2*a*d)*x)/(c*
x + d)) + 2*(4*(a*b^2*c^5 - 2*a^2*b*c^4*d + a^3*c^3*d^2)*x^3 + (8*a*b^2*c^4*d - 29*a^2*b*c^3*d^2 + 18*a^3*c^2*
d^3)*x^2 + (4*a*b^2*c^3*d^2 - 19*a^2*b*c^2*d^3 + 12*a^3*c*d^4)*x)*sqrt((a*x + b)/x))/(a^2*b^2*c^6*d^2 - 2*a^3*
b*c^5*d^3 + a^4*c^4*d^4 + (a^2*b^2*c^8 - 2*a^3*b*c^7*d + a^4*c^6*d^2)*x^2 + 2*(a^2*b^2*c^7*d - 2*a^3*b*c^6*d^2
+ a^4*c^5*d^3)*x), -1/4*((35*a^2*b^2*c^2*d^3 - 56*a^3*b*c*d^4 + 24*a^4*d^5 + (35*a^2*b^2*c^4*d - 56*a^3*b*c^3
*d^2 + 24*a^4*c^2*d^3)*x^2 + 2*(35*a^2*b^2*c^3*d^2 - 56*a^3*b*c^2*d^3 + 24*a^4*c*d^4)*x)*sqrt(d/(b*c - a*d))*a
rctan(-(b*c - a*d)*x*sqrt(d/(b*c - a*d))*sqrt((a*x + b)/x)/(a*d*x + b*d)) - 2*(b^3*c^3*d^2 + 4*a*b^2*c^2*d^3 -
11*a^2*b*c*d^4 + 6*a^3*d^5 + (b^3*c^5 + 4*a*b^2*c^4*d - 11*a^2*b*c^3*d^2 + 6*a^3*c^2*d^3)*x^2 + 2*(b^3*c^4*d
+ 4*a*b^2*c^3*d^2 - 11*a^2*b*c^2*d^3 + 6*a^3*c*d^4)*x)*sqrt(a)*log(2*a*x - 2*sqrt(a)*x*sqrt((a*x + b)/x) + b)
- (4*(a*b^2*c^5 - 2*a^2*b*c^4*d + a^3*c^3*d^2)*x^3 + (8*a*b^2*c^4*d - 29*a^2*b*c^3*d^2 + 18*a^3*c^2*d^3)*x^2 +
(4*a*b^2*c^3*d^2 - 19*a^2*b*c^2*d^3 + 12*a^3*c*d^4)*x)*sqrt((a*x + b)/x))/(a^2*b^2*c^6*d^2 - 2*a^3*b*c^5*d^3
+ a^4*c^4*d^4 + (a^2*b^2*c^8 - 2*a^3*b*c^7*d + a^4*c^6*d^2)*x^2 + 2*(a^2*b^2*c^7*d - 2*a^3*b*c^6*d^2 + a^4*c^5
*d^3)*x), -1/4*((35*a^2*b^2*c^2*d^3 - 56*a^3*b*c*d^4 + 24*a^4*d^5 + (35*a^2*b^2*c^4*d - 56*a^3*b*c^3*d^2 + 24*
a^4*c^2*d^3)*x^2 + 2*(35*a^2*b^2*c^3*d^2 - 56*a^3*b*c^2*d^3 + 24*a^4*c*d^4)*x)*sqrt(d/(b*c - a*d))*arctan(-(b*
c - a*d)*x*sqrt(d/(b*c - a*d))*sqrt((a*x + b)/x)/(a*d*x + b*d)) - 4*(b^3*c^3*d^2 + 4*a*b^2*c^2*d^3 - 11*a^2*b*
c*d^4 + 6*a^3*d^5 + (b^3*c^5 + 4*a*b^2*c^4*d - 11*a^2*b*c^3*d^2 + 6*a^3*c^2*d^3)*x^2 + 2*(b^3*c^4*d + 4*a*b^2*
c^3*d^2 - 11*a^2*b*c^2*d^3 + 6*a^3*c*d^4)*x)*sqrt(-a)*arctan(sqrt(-a)*sqrt((a*x + b)/x)/a) - (4*(a*b^2*c^5 - 2
*a^2*b*c^4*d + a^3*c^3*d^2)*x^3 + (8*a*b^2*c^4*d - 29*a^2*b*c^3*d^2 + 18*a^3*c^2*d^3)*x^2 + (4*a*b^2*c^3*d^2 -
19*a^2*b*c^2*d^3 + 12*a^3*c*d^4)*x)*sqrt((a*x + b)/x))/(a^2*b^2*c^6*d^2 - 2*a^3*b*c^5*d^3 + a^4*c^4*d^4 + (a^
2*b^2*c^8 - 2*a^3*b*c^7*d + a^4*c^6*d^2)*x^2 + 2*(a^2*b^2*c^7*d - 2*a^3*b*c^6*d^2 + a^4*c^5*d^3)*x)]
________________________________________________________________________________________
Sympy [F]
time = 0.00, size = 0, normalized size = 0.00 \begin {gather*} \int \frac {x^{3}}{\sqrt {a + \frac {b}{x}} \left (c x + d\right )^{3}}\, dx \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c+d/x)**3/(a+b/x)**(1/2),x)
[Out]
Integral(x**3/(sqrt(a + b/x)*(c*x + d)**3), x)
________________________________________________________________________________________
Giac [B] Leaf count of result is larger than twice the leaf count of optimal. 901 vs.
\(2 (222) = 444\).
time = 2.23, size = 901, normalized size = 3.60 \begin {gather*} \frac {{\left (35 \, a^{\frac {3}{2}} b^{2} c^{2} d^{2} \arctan \left (\frac {\sqrt {a} d}{\sqrt {b c d - a d^{2}}}\right ) - 56 \, a^{\frac {5}{2}} b c d^{3} \arctan \left (\frac {\sqrt {a} d}{\sqrt {b c d - a d^{2}}}\right ) + 24 \, a^{\frac {7}{2}} d^{4} \arctan \left (\frac {\sqrt {a} d}{\sqrt {b c d - a d^{2}}}\right ) - 2 \, \sqrt {b c d - a d^{2}} b^{3} c^{3} \log \left ({\left | b \right |}\right ) - 8 \, \sqrt {b c d - a d^{2}} a b^{2} c^{2} d \log \left ({\left | b \right |}\right ) + 22 \, \sqrt {b c d - a d^{2}} a^{2} b c d^{2} \log \left ({\left | b \right |}\right ) - 12 \, \sqrt {b c d - a d^{2}} a^{3} d^{3} \log \left ({\left | b \right |}\right ) + 13 \, \sqrt {b c d - a d^{2}} a^{2} b c d^{2} - 10 \, \sqrt {b c d - a d^{2}} a^{3} d^{3}\right )} \mathrm {sgn}\left (x\right )}{4 \, {\left (\sqrt {b c d - a d^{2}} a^{\frac {3}{2}} b^{2} c^{6} - 2 \, \sqrt {b c d - a d^{2}} a^{\frac {5}{2}} b c^{5} d + \sqrt {b c d - a d^{2}} a^{\frac {7}{2}} c^{4} d^{2}\right )}} + \frac {{\left (35 \, b^{2} c^{2} d^{2} - 56 \, a b c d^{3} + 24 \, a^{2} d^{4}\right )} \arctan \left (-\frac {{\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )} c + \sqrt {a} d}{\sqrt {b c d - a d^{2}}}\right )}{4 \, {\left (b^{2} c^{6} \mathrm {sgn}\left (x\right ) - 2 \, a b c^{5} d \mathrm {sgn}\left (x\right ) + a^{2} c^{4} d^{2} \mathrm {sgn}\left (x\right )\right )} \sqrt {b c d - a d^{2}}} + \frac {13 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{3} \sqrt {a} b^{2} c^{3} d^{2} - 40 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{3} a^{\frac {3}{2}} b c^{2} d^{3} + 24 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{3} a^{\frac {5}{2}} c d^{4} + 7 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{2} a b^{2} c^{2} d^{3} - 56 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{2} a^{2} b c d^{4} + 40 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{2} a^{3} d^{5} + 11 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )} \sqrt {a} b^{3} c^{2} d^{3} - 60 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )} a^{\frac {3}{2}} b^{2} c d^{4} + 40 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )} a^{\frac {5}{2}} b d^{5} - 13 \, a b^{3} c d^{4} + 10 \, a^{2} b^{2} d^{5}}{4 \, {\left (\sqrt {a} b^{2} c^{6} \mathrm {sgn}\left (x\right ) - 2 \, a^{\frac {3}{2}} b c^{5} d \mathrm {sgn}\left (x\right ) + a^{\frac {5}{2}} c^{4} d^{2} \mathrm {sgn}\left (x\right )\right )} {\left ({\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )}^{2} c + 2 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )} \sqrt {a} d + b d\right )}^{2}} + \frac {\sqrt {a x^{2} + b x}}{a c^{3} \mathrm {sgn}\left (x\right )} + \frac {{\left (b c + 6 \, a d\right )} \log \left ({\left | 2 \, {\left (\sqrt {a} x - \sqrt {a x^{2} + b x}\right )} \sqrt {a} + b \right |}\right )}{2 \, a^{\frac {3}{2}} c^{4} \mathrm {sgn}\left (x\right )} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate(1/(c+d/x)^3/(a+b/x)^(1/2),x, algorithm="giac")
[Out]
1/4*(35*a^(3/2)*b^2*c^2*d^2*arctan(sqrt(a)*d/sqrt(b*c*d - a*d^2)) - 56*a^(5/2)*b*c*d^3*arctan(sqrt(a)*d/sqrt(b
*c*d - a*d^2)) + 24*a^(7/2)*d^4*arctan(sqrt(a)*d/sqrt(b*c*d - a*d^2)) - 2*sqrt(b*c*d - a*d^2)*b^3*c^3*log(abs(
b)) - 8*sqrt(b*c*d - a*d^2)*a*b^2*c^2*d*log(abs(b)) + 22*sqrt(b*c*d - a*d^2)*a^2*b*c*d^2*log(abs(b)) - 12*sqrt
(b*c*d - a*d^2)*a^3*d^3*log(abs(b)) + 13*sqrt(b*c*d - a*d^2)*a^2*b*c*d^2 - 10*sqrt(b*c*d - a*d^2)*a^3*d^3)*sgn
(x)/(sqrt(b*c*d - a*d^2)*a^(3/2)*b^2*c^6 - 2*sqrt(b*c*d - a*d^2)*a^(5/2)*b*c^5*d + sqrt(b*c*d - a*d^2)*a^(7/2)
*c^4*d^2) + 1/4*(35*b^2*c^2*d^2 - 56*a*b*c*d^3 + 24*a^2*d^4)*arctan(-((sqrt(a)*x - sqrt(a*x^2 + b*x))*c + sqrt
(a)*d)/sqrt(b*c*d - a*d^2))/((b^2*c^6*sgn(x) - 2*a*b*c^5*d*sgn(x) + a^2*c^4*d^2*sgn(x))*sqrt(b*c*d - a*d^2)) +
1/4*(13*(sqrt(a)*x - sqrt(a*x^2 + b*x))^3*sqrt(a)*b^2*c^3*d^2 - 40*(sqrt(a)*x - sqrt(a*x^2 + b*x))^3*a^(3/2)*
b*c^2*d^3 + 24*(sqrt(a)*x - sqrt(a*x^2 + b*x))^3*a^(5/2)*c*d^4 + 7*(sqrt(a)*x - sqrt(a*x^2 + b*x))^2*a*b^2*c^2
*d^3 - 56*(sqrt(a)*x - sqrt(a*x^2 + b*x))^2*a^2*b*c*d^4 + 40*(sqrt(a)*x - sqrt(a*x^2 + b*x))^2*a^3*d^5 + 11*(s
qrt(a)*x - sqrt(a*x^2 + b*x))*sqrt(a)*b^3*c^2*d^3 - 60*(sqrt(a)*x - sqrt(a*x^2 + b*x))*a^(3/2)*b^2*c*d^4 + 40*
(sqrt(a)*x - sqrt(a*x^2 + b*x))*a^(5/2)*b*d^5 - 13*a*b^3*c*d^4 + 10*a^2*b^2*d^5)/((sqrt(a)*b^2*c^6*sgn(x) - 2*
a^(3/2)*b*c^5*d*sgn(x) + a^(5/2)*c^4*d^2*sgn(x))*((sqrt(a)*x - sqrt(a*x^2 + b*x))^2*c + 2*(sqrt(a)*x - sqrt(a*
x^2 + b*x))*sqrt(a)*d + b*d)^2) + sqrt(a*x^2 + b*x)/(a*c^3*sgn(x)) + 1/2*(b*c + 6*a*d)*log(abs(2*(sqrt(a)*x -
sqrt(a*x^2 + b*x))*sqrt(a) + b))/(a^(3/2)*c^4*sgn(x))
________________________________________________________________________________________
Mupad [B]
time = 5.48, size = 2890, normalized size = 11.56 \begin {gather*} \text {Too large to display} \end {gather*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int(1/((a + b/x)^(1/2)*(c + d/x)^3),x)
[Out]
(log((d^3*(a*d - b*c)^5)^(1/2)*(a + b/x)^(1/2) - a^3*d^4 + b^3*c^3*d - 3*a*b^2*c^2*d^2 + 3*a^2*b*c*d^3)*(d^3*(
a*d - b*c)^5)^(1/2)*(3*a^2*d^2 + (35*b^2*c^2)/8 - 7*a*b*c*d))/(b^5*c^9 - a^5*c^4*d^5 + 5*a^4*b*c^5*d^4 + 10*a^
2*b^3*c^7*d^2 - 10*a^3*b^2*c^6*d^3 - 5*a*b^4*c^8*d) - ((b*(a + b/x)^(5/2)*(12*a^2*d^4 + 4*b^2*c^2*d^2 - 19*a*b
*c*d^3))/(4*a*c^3*(a*d - b*c)^2) - ((a + b/x)^(1/2)*(4*b^4*c^3 - 12*a^3*b*d^3 + 25*a^2*b^2*c*d^2 - 12*a*b^3*c^
2*d))/(4*a*c^3*(a*d - b*c)) + (d*(a + b/x)^(3/2)*(8*b^4*c^3 - 24*a^3*b*d^3 + 56*a^2*b^2*c*d^2 - 37*a*b^3*c^2*d
))/(4*c^3*(a^2*d - a*b*c)*(a*d - b*c)))/((a + b/x)^2*(3*a*d^2 - 2*b*c*d) - (a + b/x)*(3*a^2*d^2 + b^2*c^2 - 4*
a*b*c*d) - d^2*(a + b/x)^3 + a^3*d^2 + a*b^2*c^2 - 2*a^2*b*c*d) - (log((d^3*(a*d - b*c)^5)^(1/2)*(a + b/x)^(1/
2) + a^3*d^4 - b^3*c^3*d + 3*a*b^2*c^2*d^2 - 3*a^2*b*c*d^3)*(d^3*(a*d - b*c)^5)^(1/2)*(24*a^2*d^2 + 35*b^2*c^2
- 56*a*b*c*d))/(8*(b^5*c^9 - a^5*c^4*d^5 + 5*a^4*b*c^5*d^4 + 10*a^2*b^3*c^7*d^2 - 10*a^3*b^2*c^6*d^3 - 5*a*b^
4*c^8*d)) - (atan((((((a + b/x)^(1/2)*(1152*a^6*b^2*d^9 + 16*b^8*c^6*d^3 + 128*a*b^7*c^5*d^4 - 4800*a^5*b^3*c*
d^8 + 1129*a^2*b^6*c^4*d^5 - 5136*a^3*b^5*c^3*d^6 + 7520*a^4*b^4*c^2*d^7))/(8*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4*
a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2)) - (((4*a*b^8*c^13*d^2 + 4*a^2*b^7*c^12*d^3 - 45*a^3*b^6*
c^11*d^4 + 74*a^4*b^5*c^10*d^5 - 49*a^5*b^4*c^9*d^6 + 12*a^6*b^3*c^8*d^7)/(a^2*b^4*c^13 + a^6*c^9*d^4 - 4*a^3*
b^3*c^12*d - 4*a^5*b*c^10*d^3 + 6*a^4*b^2*c^11*d^2) - ((a + b/x)^(1/2)*(6*a*d + b*c)*(64*a^2*b^7*c^13*d^2 - 38
4*a^3*b^6*c^12*d^3 + 896*a^4*b^5*c^11*d^4 - 1024*a^5*b^4*c^10*d^5 + 576*a^6*b^3*c^9*d^6 - 128*a^7*b^2*c^8*d^7)
)/(16*c^4*(a^3)^(1/2)*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4*a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2)))*(
6*a*d + b*c))/(2*c^4*(a^3)^(1/2)))*(6*a*d + b*c)*1i)/(2*c^4*(a^3)^(1/2)) + ((((a + b/x)^(1/2)*(1152*a^6*b^2*d^
9 + 16*b^8*c^6*d^3 + 128*a*b^7*c^5*d^4 - 4800*a^5*b^3*c*d^8 + 1129*a^2*b^6*c^4*d^5 - 5136*a^3*b^5*c^3*d^6 + 75
20*a^4*b^4*c^2*d^7))/(8*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4*a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2))
+ (((4*a*b^8*c^13*d^2 + 4*a^2*b^7*c^12*d^3 - 45*a^3*b^6*c^11*d^4 + 74*a^4*b^5*c^10*d^5 - 49*a^5*b^4*c^9*d^6 +
12*a^6*b^3*c^8*d^7)/(a^2*b^4*c^13 + a^6*c^9*d^4 - 4*a^3*b^3*c^12*d - 4*a^5*b*c^10*d^3 + 6*a^4*b^2*c^11*d^2) +
((a + b/x)^(1/2)*(6*a*d + b*c)*(64*a^2*b^7*c^13*d^2 - 384*a^3*b^6*c^12*d^3 + 896*a^4*b^5*c^11*d^4 - 1024*a^5*b
^4*c^10*d^5 + 576*a^6*b^3*c^9*d^6 - 128*a^7*b^2*c^8*d^7))/(16*c^4*(a^3)^(1/2)*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4*
a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2)))*(6*a*d + b*c))/(2*c^4*(a^3)^(1/2)))*(6*a*d + b*c)*1i)/(
2*c^4*(a^3)^(1/2)))/((216*a^5*b^3*d^9 + (35*b^8*c^5*d^4)/2 - (49*a*b^7*c^4*d^5)/8 - 810*a^4*b^4*c*d^8 - (1877*
a^2*b^6*c^3*d^6)/4 + 1044*a^3*b^5*c^2*d^7)/(a^2*b^4*c^13 + a^6*c^9*d^4 - 4*a^3*b^3*c^12*d - 4*a^5*b*c^10*d^3 +
6*a^4*b^2*c^11*d^2) + ((((a + b/x)^(1/2)*(1152*a^6*b^2*d^9 + 16*b^8*c^6*d^3 + 128*a*b^7*c^5*d^4 - 4800*a^5*b^
3*c*d^8 + 1129*a^2*b^6*c^4*d^5 - 5136*a^3*b^5*c^3*d^6 + 7520*a^4*b^4*c^2*d^7))/(8*(a^2*b^4*c^10 + a^6*c^6*d^4
- 4*a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2)) - (((4*a*b^8*c^13*d^2 + 4*a^2*b^7*c^12*d^3 - 45*a^3*
b^6*c^11*d^4 + 74*a^4*b^5*c^10*d^5 - 49*a^5*b^4*c^9*d^6 + 12*a^6*b^3*c^8*d^7)/(a^2*b^4*c^13 + a^6*c^9*d^4 - 4*
a^3*b^3*c^12*d - 4*a^5*b*c^10*d^3 + 6*a^4*b^2*c^11*d^2) - ((a + b/x)^(1/2)*(6*a*d + b*c)*(64*a^2*b^7*c^13*d^2
- 384*a^3*b^6*c^12*d^3 + 896*a^4*b^5*c^11*d^4 - 1024*a^5*b^4*c^10*d^5 + 576*a^6*b^3*c^9*d^6 - 128*a^7*b^2*c^8*
d^7))/(16*c^4*(a^3)^(1/2)*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4*a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2)
))*(6*a*d + b*c))/(2*c^4*(a^3)^(1/2)))*(6*a*d + b*c))/(2*c^4*(a^3)^(1/2)) - ((((a + b/x)^(1/2)*(1152*a^6*b^2*d
^9 + 16*b^8*c^6*d^3 + 128*a*b^7*c^5*d^4 - 4800*a^5*b^3*c*d^8 + 1129*a^2*b^6*c^4*d^5 - 5136*a^3*b^5*c^3*d^6 + 7
520*a^4*b^4*c^2*d^7))/(8*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4*a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2))
+ (((4*a*b^8*c^13*d^2 + 4*a^2*b^7*c^12*d^3 - 45*a^3*b^6*c^11*d^4 + 74*a^4*b^5*c^10*d^5 - 49*a^5*b^4*c^9*d^6 +
12*a^6*b^3*c^8*d^7)/(a^2*b^4*c^13 + a^6*c^9*d^4 - 4*a^3*b^3*c^12*d - 4*a^5*b*c^10*d^3 + 6*a^4*b^2*c^11*d^2) +
((a + b/x)^(1/2)*(6*a*d + b*c)*(64*a^2*b^7*c^13*d^2 - 384*a^3*b^6*c^12*d^3 + 896*a^4*b^5*c^11*d^4 - 1024*a^5*
b^4*c^10*d^5 + 576*a^6*b^3*c^9*d^6 - 128*a^7*b^2*c^8*d^7))/(16*c^4*(a^3)^(1/2)*(a^2*b^4*c^10 + a^6*c^6*d^4 - 4
*a^3*b^3*c^9*d - 4*a^5*b*c^7*d^3 + 6*a^4*b^2*c^8*d^2)))*(6*a*d + b*c))/(2*c^4*(a^3)^(1/2)))*(6*a*d + b*c))/(2*
c^4*(a^3)^(1/2))))*(6*a*d + b*c)*1i)/(c^4*(a^3)^(1/2))
________________________________________________________________________________________